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Abstract 

A scalar field, ?7(r), is coupled to the skyrmion. Its classical value in the vacuum ("condensate") 
reduces to the pion decay constant, f n , being thus proportional to the chiral condensate < qq >. 
A quadratic coupling of the scalar field to the pion kinetic Lagrangian term, reminiscent from the 
linear sigma model, is considered. Its mass is an additional free parameter in the model associated 
to the lightest scalar isoscalar meson. Solutions for the two resulting coupled differential equations 
are found in several cases. The nucleon, as a topological soliton, "digs a hole" in the vacuum 
making the value of the scalar condensate to be close to zero inside the skyrmion. Some of the 
observables, as for example the baryon masses, may have their values closer to the experimental 
ones. Derrick's stability analysis is applied to the model. 
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1 Introduction 



The description of processes involving strong interactions (mainly at low and intermediate energies) 
is very difficult in the frame of the Quantum Chromodynamics (QCD) due to its non abelian color 
and flavor structures and strong coupling constants. Thus effective models, like the Skyrme chiral 
soliton model (Skyrmion model), are constructed in such a way as to respect general properties from 
the more fundamental theory, such as chiral symmetry and its spontaneous breakdown p]. 

The topological soliton model was formulated by Skyrme decades ago [^j possessing a conserved 
topological charge identified with the baryonic number. A deeper discussion about the real existence 
of minimum energy configurations in this model is found for example in (3|. The model encompasses 
chiral symmetry and its spontaneous breakdown and has further deeper motivation based on the large 
N c QCD analysis |3|, ^]. The quantization of the rotational collective coordinates endows the soliton 
with spin half. It strengthens the interpretation of this soliton as a nucleon as long as it is possible 
to obtain predictions for hadronic properties @. An introductory study of the inclusion of quantum 
and thermal effects was done in [JjJ. 

Chiral spontaneous symmetry breaking is believed to be realized in the Nambu-Goldstone mode: 
the pion decay constant, /„-, is different from zero and the pion mass is small in the hadronic scale. 
f n can be seen as a "sigma condensate" of the linear sigma model. This sigma classical value is 
also identified to the chiral radius which is the constraint for the non linear realization of chiral 
symmetry as well. The non zero pion mass is usually assumed to break the chiral invariance of the 
Lagrangian explicitly, although only slightly. The vacuum is expected to be fulfilled with the scalar 
quark condensate (< qq >) which is proportional to the pion decay constant. To the lowest order, the 
Gell-Mann-Oakes-Renner relation, which relates QCD to hadronic degrees of freedom, reads: 



The experimental value of the pion decay constant is nearly 93 MeV. However, in the context of 
the Skyrmion model, it may reach very low values (such as 54 MeV) as long as one takes as input 
the physical values of nucleon and delta masses. This may suggest by itself that f n can acquire 
different values inside and outside the nucleon. The constant value of the chiral radius is one trivial 
(homogeneous) solution of the linear sigma model. We show in this work that there are other non 
trivial (non homogeneous) solutions in the presence of a nucleon, i.e., it occurs the formation of a 
"hole" in the vacuum due to the presence of the nucleon in the frame of the Skyrmion model. 



2 < qq > m q 
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The coupling to QCD scalar degrees of freedom has already been envisaged in the frame of the 
Skyrmion model but in different approaches from that exploited here. Schechter and collaborators || 
and afterwards Waindzoch and Wambach § investigated such kind of effect by coupling a scalar field, 
representing a gluon degree of freedom. The value obtained for that "gluon condensate" inside the 
nucleon was much higher than in the vacuum, thus forming a kind of "bag". This strengthens some 
qualitative theoretical basis of the so called Bag Models [HI. In particular, for light glueballs, the 
bag was found to be very shallow (8|. These results are in agreement with the stabilizing mechanism 
discussed in . However the QCD vacuum is very complex and it is expected to contain, as discussed 
above, the quark-antiquark scalar condensate which is a manifestation of the spontaneously broken 
chiral symmetry. Thus a more detailed investigation of the contribution of this complex system for 
the structure of the nucleon is needed as well as the opposite, i.e., the effect of the nucleon on the 



vacuum. A more complex system was considered in [12] where a scalar meson field was also included. 
We have earlier considered the chiral radius in the presence of hadrons as a dynamical variable 



13, 14]. In this work we consider a dynamical scalar degree of freedom associated with the scalar 



chiral condensate. This is effectively implemented by the coupling of the Skyrmion to a self interacting 
scalar field. As shown is sections II and III, the scalar condensate ("sigma condensate") is constrained 
by the form of the potential which endows the model with the spontaneous symmetry breaking: in 
the vacuum it assumes the usual value of the pion decay constant, rj vacuum = 93 MeV. In the 
next sections we show the Lagrangian density and equations of the fields of this extended Skyrmion 
model. There are two coupled differential equations: for the soliton profile function, F(r), and for 
the scalar isoscalar field, r](r). The numerical solutions are investigated for some cases and hadronic 
observables are calculated. We investigate the effect of the magnitude of the mass of this scalar field 
- the additional free parameter of our model - privileging values between 560 MeV and 980 MeV, 



which correspond respectively to the sigma (eventually the chiral partner of the pion) (T^, 16, 17 
and to the /o(980) meson fL8fl . The latter is one of the lightest well known scalar isoscalar mesons. 
The analysis of stability of the equations by means of dimensional arguments is also discussed in the 
section IV (Derrick's analysis). In section V we present our final comments, also discussing some ways 
of improving and extending the model. 
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2 Scalar field coupling 

We do not describe the Skyrmion model at length since this is done extensively in the literature. 
One natural way of imposing a coupling of a scalar field (such as the chiral partner of the pion) to 
the skyrmion model is to consider the linear realization of the sigma model. The degree of freedom 
of the sigma may be replaced by a variable (dynamical) chiral radius, r/(r), which becomes the new 
dynamical degree of freedom: 



fir^V = Vn 2 + a 2 = v + e, 



where v = f n is the classical value of the scalar field and e its fluctuation [19]. However, the resulting 
coupled model (roughly the linear sigma model) has no stable solitonic solution with non trivial 
topology as discussed later in this paper. Thus we include, as it is usually done in the standard 
Skyrmion model, the fourth order stabilizing term for the pionic sector. The modified Skyrmion 
Lagrangian density reads: 

C = ^^Tr(d^d^U) + ^Tr ([rfdpU, tfd v uf} + ^ 2 Uv(r) (TrU + Trrf - 2) + C v , (2) 

where the first term is the Lagrangian density of the linear sigma model and the second one is 
introduced in order to stabilize the soliton Q. While the first term becomes scale invariant with 
the quadratic coupling, the second one is scale invariant by itself and thus no such coupling has been 
considered in this work for this term 0. fj 2 is the pion mass whose Lagrangian term breaks slightly 
the chiral symmetry and scale invariance as well, at least in this form. U is a chiral invariant matrix 
parametrized in terms of the profile function F(r) within the Hedgehog ansatz (given by 7r = r): 

U(r) = cos(F(r)) + it ■ 7r sin(F(r)). 

Finally, the scalar field self interacting Lagrangian is chosen such that one recovers a classical value 
in the vacuum due to the (chiral) spontaneous symmetry breaking: 

A = \d^{v)d^(v) - \my(r) - \\rj\r). (3) 

The scalar field mass term breaks scale invariance but it is necessary for the consistency of the model. 



The bare scalar field mass parameter is m 2 being related to the physical scalar field mass (m 2 ) by 



v 

the following relation: 



m l = \^ 2 -m 2 t ). (4) 
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As a matter of fact, it is possible to introduce a scale invariant coupling in such term. 
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The value of the condensate in the vacuum is found by the minimization of its potential. As stated 
above we consider two different possibilities for the physical particle associated with the scalar field. 
The T) field is expected to correspond to the lightest scalar-isoscalar meson. It may be a very broad 
scalar resonance in the s-wave of low energy ir — ir scattering, with a mass in the range 500 — 800 MeV 
that seems to have been observed in other processes [f|| H, 0. In this case it would be associated 
to the chiral partner of the pion eventually corresponding to a Higgs field for Hadronic interactions i. 
We consider two different values for the scalar-isoscalar particle mass in order to investigate its effects 
on nucleon properties. The coupling constant e may be related to the pion-pion scattering coupling 



as done in 21], i.e., to a particular choice of the low energy constants (LECs) of some terms from 



the Chiral Perturbation Theory expansion [23]. In the works just quoted above the authors obtained: 
e ~ 4 — 7. Large N c analysis with extrapolation to N c = 3 using sum rules of [24] leads to the 
approximated values: e ~ (7.6 — 12) jyj N c . These values are, in principle, independent of /„-. 

The coupling to the quadratic Lagrangian term of the Skyrmion model, as described above, is 
uniquely defined by the linear implementation of chiral symmetry. It is scale invariant. It is consis- 
tent with the consideration that the chiral radius is directly and unambiguosly related to the quark- 
antiquark condensate of QCD and, besides, to the pion decay constant, at least in the vacuum and at 
the tree level. For the coupling of the quartic stabilizing Skyrme term to the scalar field we found no 
scale invariant form. 

Along this work we basically assume the chiral radius to be a dynamical variable, being equivalent 
to the sigma field. In this case one is led to the following replacement for the quadratic Skyrmion 
Lagrangian term: 

U(r) -» a(r)C/(r), with a(r) = (5) 

However, this replacement in the fourth order Skyrme term yields scale non- invariant Lagrangian 
terms and unstable solutions in the physical region for the parameters of the model. The lack of 
stability was also checked by the Derrick's argument |^5| discussed in section IV. 

Finally, one can notice that the fourth order Skyrme term may be derived from the coupling of 



the sigma model with the vector meson, p M , as a hidden or massive gauge boson [22, Pq |. The above 



prescription for the function U(r) may yield a r/(r)-dependent p mass and e coupling such that scale 
invariance is mantained. This will be developed elsewhere. Although this is a very interesting idea 



Its origin, on the other hand, lies on the chiral spontaneous symmetry breaking of QCD vacuum. 
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by itself, after the elimination of the p degrees of freedom, the prescription seems to lead to other 
derivative and complicated terms which also may break scale invariance. At last it is interesting to 
preserve scale invariance in spite of the massive terms. 



3 Calculations 

In this section we investigate the solitonic solutions relative to equations of the fields of the model and 
afterwards their consequences for static nucleon observables. Adopting the usual hedgehog ansatz, 
7r = f , the modified Lagrangian is written as: 

f , 2 f 1 2 n /2 sin 2 (F) . 2 sin 2 (F) sin 4 (F) 2 , 1 
J drr \-\ F ~ ~ —pt 1 ~ + ^Ucos(F) + £,| , (6) 

where F = F(r), r\ = r/(r) and ' is the derivative with relation to the radial coordinate r. 

The Euler-Lagrange equations for this system with the change of variable r — > rf^ are given by: 



2 2sin 2 (F)\ n , . sin(2F) 2 . , nri . 4sin 3 (F)cos(F) . . 

2 + ^ + 2rrj 2 F' + S^F' 2 - r] 2 stn(2F) ^ ' v 1 {+} 

e z J e z 2e z r z 



— fi 2 r]sin(F)r 2 f n + 2rr]r]'F' = 0, 

t?" + V - r? (V + 2 "^ 2 2 (F) + m 2 ^) + p 2 cos{F)U - W = 0- (8) 

In the spatially homogeneous case of r](r) = f n (constant), the usual skyrmion equation and solutions 
are obtained for the first equation. That value for the sigma "condensate" (classical field) naturally 
emerges as a boundary solution in the limit r — > oo, i.e., in the vacuum, far from the topological 
soliton for the solutions of interest. In the limit of m 2 — ► oo the equations become uncoupled, what is 
verified numerically. 

For the scalar field equation we consider the following conditions: 

V (r - oo) = U, (9) 
rj (r — > oo) = rj'(r = 0) = 0. 

In this limit, in the vacuum, a constraint between the scalar field parameters is obtained. By mini- 
mizing the (homogeneous) scalar field potential from equation (j8|) we obtain: 

A = (10) 

f4 
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The corresponding solutions for the spatially non homogeneous case are presented later. Spontaneous 
chiral symmetry breaking occurs as long as 

A < °> (ii) 

which provides a condition for the existence of the "sigma condensate" at the classical level. Quantum 
fluctuations may considerably change this |27], They will be investigated elsewhere. 

For the equation of the profile function F(r) the usual boundary conditions still hold, i.e., winding 
number n = 1, for which one considers F(r = 0) = ir and F(r — ► oo) = ||. From expression (^|) one 
sees that it is possible to define an effective mass parameter for the field rj(r). Namely: 

</=ro?+F « + «. (12) 

Far from the nucleon this effective mass is just equal to m^, but inside the skyrmion and close to it 
the contributions from the two other terms may be relevant. 

3.1 Numerical results 

We are thus faced with two differential equations with two point-type boundary conditions. They 
are solved in an iterative way by the "shooting" method using fourth order Runge Kutta with the 
boundary conditions shown above. 

The numerical solutions for the skyrmion profile function are presented in figure 1 and, in figure 
2, the solutions for the condensate ij(r). Three different cases are shown considering f n = 93 MeV. 
The following values for the other input parameter were used for each solution: 

(f) m € = 980MeV and e = 4, 

(2) m e = 980MeV and e = 5, (13) 

(3) m £ = 560MeV and e = 4. 

From figure 1 we notice that the skyrmion radius becomes a little smaller and that inside the topological 
soliton the scalar field tends to decrease, i.e., the nucleon "digs a hole" in the vacuum otherwise fulfilled 
with the condensate. For smaller (bigger) values of the scalar field mass the "hole" - which is in fact 
akin to the "bag" of other models - becomes deeper and wider (shallower and narrower). This is a 
mechanism for bag formation equal to those studied in H ^, 12]. As a matter of fact, the Skyrmion 



model and the constituent quark model were shown to yield equivalent results in the large N c limit 



[29 1, thus the bag structure should be expected in both models. 



In fact, we find that inside the topolotical soliton the effective mass given in expression (|T^) 
should be smaller in modulus, being even positive, thus violating condition ([ll]) for the occurrence of 
spontaneous symmetry breaking. It can be seen as the restoration of the chiral symmetry inside the 
soliton (nucleon) for not high values of ml. Therefore there is a sort of a "critical value" for the scalar 
particle mass above which the chiral symmetry is restored inside the nucleon. This value depends on 
the other parameters of the model and on the level of approximation for the quantum fluctuations of 
the fields. 

Although it is not clear whether it is physically meaningful we describe the following curious 
issue. We have also investigated the possibility that at the origin, i.e., inside the nucleon, the classical 
scalar field is larger than in the vacuum. For a very particular fine tuned value of r](r —►())- which 
depends on e, m 2 and rj(r — ► 0) - we found a solution which goes asymptoticaly to /„-, as required and 
expected. We have obtained this only for large mass values. However this solution exhibits a sort of 
small amplitude oscillation, which decreases progressively with r, having large wavelength around the 
region of the topological soliton. It would correspond to a "stronger chiral symmetry breaking" inside 
the soliton. We leave the investigation of such solution for another work |3Cj. 

Strong and electromagnetic properties usually obtained for the skyrmion were calculated con- 
sidering the above scalar field coupling. The energy of the lightest baryons (nucleons) can be written 
as: 



M = 4vr jT drr 2 K (V 2 + ^^-) + (1 - cos{F))^U + ^ + ^ 2 + 

1 f sin 2 (F)F' 2 sin 4 (F)\ 1 1 
+ ^2 ~2 + 2r^ } + X~ | ' 



(14) 



where I is the moment of inertia j|] which now depends on r)(r) and E vac is the vacuum energy due 
to the dynamical chiral symmetry breaking to be subtracted: 

In table 1 we compare the values of observables obtained for the usual Skyrmion model, our model 
(for the three cases shown in figures 1 and 2) and experimental data. A higher value for e (e = 5) yield 
better values for the mass spectrum as it is well known. However almost all the other observables 
(mainly the electromagnetic properties) still present values in disagreement with the experimental 
ones. On the other hand, considering smaller values of the scalar field mass, m e = 560 MeV, a 
better agreement with the experimental values for most observables is obtained. The most remarkable 



exception is the behaviour of the axial coupling constant qa- The qa coupling is calculated in two 
ways: 1) in the limit of zero momentum of the nucleon axial form factor and 2) by means of the 
Goldberger-Treiman relation || in which case all the parameters can be calculated from the resulting 
Skyrmion solutions. For gA and g n NN a smaller value of e would be desirable, what would, in principle, 
make the values of the other observables to be in disagreement with experimental ones. 



4 Derrick's analysis of stability 

We have also performed the Derrick's analysis in order to verify the stability of the "finite sized" 
solutions of the differential equations |2f|. It consists in checking, by dimensional arguments (scale 
transformations), whether there is a scale in which the total energy of the system is finite yielding 
stable finite size solutions for the corresponding equations. One has to consider scale transformations 
for the spatial coordinate variable 

r — ► r 



J r 



a 



and for the dimensionful field 



rj(r) — > (3r)[r r ). 



In fact it is more suitable to do it for the shifted equivalent field: 9(r) = f n — rj(r). Under the above 
transformations we verify whether the energy density has a stable minimum with relation to the scale 
transformation parameters, i.e.: 

8M 5M 

^=°- w- - (16) 
fo^ >0 - w >0 - (17) 

This analysis was performed for the following cases: 

(i) Only the sigma model with the scalar field which has self interacting potential, i.e., no fourth order 
term for the pion field in the Lagrangian density (Q); 

(ii) The inclusion of the usual quartic Skyrme term exactly as it stands in expression (0); 

(iii) The coupling of the scalar field in the fourth order Skyrme term, breaking the scale invariance, 
with the ansatz of expression (|5|). This modification yields, respectively, the following modified and 
additional Lagrangian density terms: 

1 f ,,, /s 2 F' 2 s 4 \ rfWl , > 

tl = W {" ^ — + 2H ) + -^ ' <18) 



In the first case (i) there are no stable solutions for the skyrmion profile function (5 2 M/5a 2 < 0). 
It becomes clear once again the importance of the higher order derivative terms in order to keep the 
soliton against collapse, in spite of the problems with the convergence of the derivative series [[II]] . 



Case (ii) is stable provided that the chiral symmetry breaking term is small (as always assumed) 
compared to the self interacting scalar field term as it can be seen from the second equation of ( P~9| ) 
below. The conditions for stability are given by: 
5 2 M f°° . 1 



> — > / r 2 dr- 



5a 2 ' Jo ef n 

5 2 M 



>0, 



5(3 2 



> -» J r 2 dr (\ri A + //r]U(cos(F) - 1)) > 0. 



(19) 

The first condition does not depend on e neither on f n . However the second one depends on the value 
of A. This is found in numerical calculations where for very small values of A (and therefore of m^) 
the solution for ij(r) becomes so stiffer that it makes the numerical calculations more unreliable and 
the solutions eventually inexistent for the boundary conditions we consider. 



For case (iii) we have the following stability conditions: 



da 2 Jo \ f 2 \ r 2 r 4 I 4 1 r 2 

5 2 M f°° 2 f , i a ( sin 2 (F)(F') 2 sin 4 (F)\ 2 „ . '.\ 
jpT>0^ J o r^rhA^ + iee 2 ^ 4 ! K J K > + ^ j - /j, f w rj(l - cos(F)) j > 0. 

(20) 

The first of these equations depends on a scale (f n or e) for the stability and the second one depends 
on the relation between two parameters, A and e. This can be seen by means of an usual change of 
variable such as r — ► ef n r. 

With this analysis we see that the coupling of the scalar field to the fourth order Skyrme term as 
in prescription (||) is not suitable for the stability of the Skyrmion solution as it was argued before. 

5 Summary and Conclusions 

We have investigated the effect of coupling a classical scalar isoscalar field to the usual Skyrmion 
model as a degree of freedom whose corresponding quantum would be the chiral partner of the pion. 



The dependence of the observables on the values of the scalar field mass and on the fourth order 
coupling term e were also considered. Smaller values of the a mass favor better agreement for most 
of the nucleon static properties, except qa- Besides the expected solution of small values of the scalar 
condensate inside the nucleon, we found the possibility of a higher value for the classical field (and 
also for its mass) inside the nucleon, but for restricted values of the free parameters. This would 
be equivalent to a further, stronger, chiral symmetry breaking [28, |3^ ]. This feature is not yet well 



understood and is being investigated |30 |. 

This non trivial solution for the chiral radius (equivalently, for the sigma classical field) may 
suggest that the chiral symmetry is realized in the Wigner-Weyl mode inside the soliton instead of the 
Nambu-Goldstone mode, which is expected to hold in the vacuum. The scalar field mass is another 
free parameter that determines the form of the potential which presents (or not) the spontaneous 
symmetry breaking. The coupling to the nucleon allows the definition of an effective mass parameter, 
given by expression (|i~2|) , which may be even positive inside the nucleon for the solutions shown in the 
figures while outside (in the vacuum) it becomes negative yielding the "mexican hat" -type potential 
for T). 

In this work the pion mass term was considered to break explicitely chiral symmetry and to have 
the more standard form. However this form is not fully determined unless for the requirement of 
yielding the leading terms of the axial current divergence - PC AC. Furthermore its coupling to the 
scalar field is not at all uniquely determined and it may be expected to be related to the trace anomaly 
of QCD [^] and certainly to the chiral SSB. Alternative schemes for that will be shown elsewhere. 

As discussed above there is a procedure for obtaining the fourth order stabilizing term from a 
coupling to the p meson field (and also to its axial chiral partner A\). It is based on the elimination 
of the vector meson degrees of freedom by performing the infinite mass limit and keeping the KSRF 



relation [26|: m p = 2g p ^ T1T f^. However, we have consistently considered that inside the nucleon f w 
becomes a dynamical degree of freedom whose value may be seemingly reduced. Therefore the limit 
of large rho mass is not necessarily reliable anymore, thus imposing difficulties in order to relate one 
picture (with the p) to another (without it). This will be shown in detail elsewhere. 
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Figure captions 
Figure 1. 

Profile function F(r) for uncoupled skyrmion (solid line) and taking into account the coupling to 
the condensate with input parameters: e=4 and m t = 980 MeV (dashed line), e=5 and m e = 980 MeV 
(dotted- dashed line) and e=4 with m t = 560 MeV (dotted line). 

Figure 2. 

Condensate rj(r) in the presence of the soliton for the three cases with coupling of figure 1. 
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Table 1: Values of the observables calculated for the free skyrmion model and the other three cases 
discussed in the text (with = 93 MeV). g A is calculated in two ways, from the zero momentum 
nucleon axial form factor as well as from the Goldberger-Treiman relation. 



observable 


free 


e=4, 


e=5, 


e=4, 


experiments 




skyrmion 


m e = 980MeV 


m e = 980MeV 


m e = 560 MeV 




M N (MeV) 


1819 


1578 


1282 


1436 


939 


M A (MeV) 


2069 


1876 


1830 


1722 


1210 


A (MeV" 1 ) 


166.7 


198 


366 


191 


— 


ci-term 


49.6 


62 


51 


54 


50 ± 20 


i 

< 4 >I(fm) 


.51 


.51 


.43 


.56 


.72 


i 

< r| >£(fm) 


.87 


.80 


.71 


.82 


.88 


i 

< r M >B ( fm ) 


.82 


.53 


.44 


.56 


.56 


< r 2 M >l (fm) 


.80 


.62 


.55 


.67 


.64 


< r| >| (fm) 


.86 


.67 


.59 


.71 


.70 


i 

< r\ >l (fm) 


-.24 


-.19 


-.16 


-.18 


-.12 


i 

< r 2 M >| (fm) 


.66 


.61 


.53 


.63 


.86±.06 


i 

< r 2 M >ri (fm) 


.68 


.64 


.59 


.65 


.88±.07 


H B (magn.) 


.34 


.35 


.37 


.36 


.44 


[i v (magn.) 


3.64 


2.66 


1.17 


2.5 


2.35 


[x p (magn.) 


3.97 


3.01 


1.54 


2.9 


2.79 


fi n (magn.) 


-3.30 


-2.31 


-.80 


-2.1 


-1.91 


HA++(magn.) 


7.55 


5.84 


3.21 


5.6 


4.7- 6.7 


H NA (magn.) 


2.66 


2.24 


1.21 


2.32 


3.3 


g A (g A GT rel.) 


1.01(1.01) 


.76(0.84) 


.49(.54) 


.78(.83) 


1.23 


gir-rrN 


20.3 


14.5 


7.5 


12.7 


13.5 


< r\ >2 fm 


.49 


.86 


.85 


.48 


.68 


< r 2 . >^fm 


.6 


.50 


.37 


.28 


.54 



15 



